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Summary 


Although stringers are used primarily as stiffen- 
ers, they also can make damaged structures fail-safe 
or damage tolerant. Assessment of the damage toler- 
ance of structures weakened by cracks is aided by 
knowledge of stress-intensity factors. In this pa- 
per, the stress-intensity factors are determined for 
a cracked orthotropic sheet adhesively bonded to an 
orthotropic stringer where the adhesive layer is al- 
lowed to debond in accordance with an adhesive fail- 
ure criterion. It is assumed that the debond penetra- 
tion under the stringer is small compared with the 
width of the stringer. Based on this assumption, the 
stringer is modeled as a semi-infinite sheet bonded to 
an infinite sheet with a finite through-the-thickness 
crack in the infinite sheet. Modeling the adhesive 
layer with a nonlinear elastic stress-strain curve is 
also considered. Since the stringer is modeled as a 
semi-infinite sheet, the solution is most appropriate 
for a crack tip located near a stringer edge. Both the 
sheet and the stringer are treated as homogeneous, 
orthotropic media which are representative of many 
fiber-reinforced composite materials. It is assumed 
that the sheet and the stringer are in a state of plane 
stress and the adhesive is in pure shear. By using 
Green’s functions and the complex variable theory of 
orthotropic elasticity developed by Lekhnitskii, a set 
of integral equations is obtained. The integral equa- 
tions are replaced by an equivalent set of algebraic 
equations, which are solved to obtain the shear-stress 
distribution in the adhesive layer. Using these adhe- 
sive stresses, the crack-tip stress- intensity factors are 
found. 

The effect of adhesive failure on the stress- 
intensity factors is examined. Debonding of the ad- 
hesive layer is found to cause an increase in the stress- 
intensity factors compared with the solution for an 
intact bond. The larger the debond area, the less ef- 
fective the stringer is in reducing the stress-intensity 
factors. A debond area with a width and height 
greater than or equal to the half-crack length vir- 
tually eliminates the effect of the stringer. 

To predict the debond areas, a rupture criterion 
based on the combined adhesive stresses is intro- 
duced. When the crack is not under the stringer, the 
majority of the debond growth is found to be along 
the edge of the stringer. When the crack tip is be- 
neath the stringer, the debond first grows to the end 
of the crack; the majority of the subsequent growth 
is along the edge of the stringer. The stress level re- 
quired to start debonding in the adhesive decreases as 
the crack tip is moved further beneath the stringer. 


The effect of debonding and a nonlinear stress- 
strain curve for the adhesive are also examined. With 
the nonlinear adhesive, the debond initiates at higher 
applied stress levels than with the linear adhesive. 
The debond areas are slightly smaller than those 
found for the linear adhesive cases. 

Debonding with the linear elastic adhesive causes 
the largest increase in the stress-intensity factors 
when compared with the intact bond case, where the 
larger the debond area, the greater the increase in 
the stress-intensity factors. Compared with the lin- 
ear adhesive solution, modeling a nonlinear adhesive 
causes the stress-intensity factor to increase when 
the adhesive bond is assumed to remain intact, but 
causes the stress-intensity factor to decrease when 
debonding is included in the model. 

Introduction 

Fibrous composite materials are light, stiff, and 
strong; therefore, they have great potential for re- 
ducing weight in aircraft structures. In general, the 
structural configurations used in composite airplanes 
have been similar to the sheet-stringer construction 
currently used in metal airplanes. In metal airplanes, 
stringers are effective in making damaged structures 
fail-safe or damage tolerant. Because of this wide 
use, the interaction of a through-the-thickness crack 
and a stringer is an important problem and has been 
investigated by many authors. 

A brief review of some of the work dealing with 
the problem of a stringer bonded to a sheet con- 
taining a crack begins with the 1959 investigation 
of Sanders (ref. 1), in which he treated the problem 
of a crack positioned symmetrically beneath a con- 
tinuously attached stringer. He employed a shear- 
lag line stringer and assumed the sheet to be in- 
extensible in the direction parallel to the crack to 
obtain the stringer stress-concentration factor and 
the crack-tip stress-intensity factor. In 1965, Greif 
and Sanders (ref. 2) removed Sanders’ previous inex- 
tensionality assumption and symmetry requirements. 
They concluded that the stringer-induced reduction 
in the stress-intensity factor was a localized effect 
and that the stringer stress-concentration factor was 
largely insensitive to sheet stiffness in the direction 
parallel to the crack. 

In 1974, Arin (ref. 3) examined the effect of a 
partially debonded, infinite stringer on the stress- 
intensity factor at the crack tip. He assumed that 
the stringer was adhesively bonded to an isotropic 
sheet along a line perpendicular to the crack, and 
he concluded that the stringer was ineffective unless 
it was quite close to the crack tip and the length 
of the debond was less than twice the crack length. 
However, since Arin was modeling the stringer as a 



line, he was not able to account for partial debond- 
ing across the width of the stringer. Anderson, Hsu, 
and McGee (ref. 4) considered the growth of a crack 
initiating at a hole in an isotropic sheet reinforced 
by a bonded doubler in 1975. Their results, based on 
a two-dimensional finite-element analysis, indicated 
that the stress-intensity factor was unusually insensi- 
tive to crack length and that adhesive debond zones 
tended to follow the tip of the crack with little sub- 
sequent growth perpendicular to the crack line. 

In 1978, Anderson, Chu, and McGee (ref. 5) as- 
sessed the growth characteristics of a fatigue crack 
approaching and growing beneath an adhesively 
bonded doubler. In this work, a two-dimensional 
finite-element analysis was used to compute the 
stress-intensity factor as a function of crack length 
for linear and nonlinear representations of the adhe- 
sive. The nonlinear representation of the adhesive 
predicted debonded areas that agreed extremely well 
with the experimentally observed debond in a bonded 
joint with aluminum adherends. The nonlinear rep- 
resentation of the adhesive was also found to reduce 
the stress-intensity factor relative to the solution for 
the linear elastic adhesive when debonding was in- 
cluded in both models of the adhesive layer. 

In 1981, using a complex-variable approach, 
Norris (ref. 6) represented the finite-width stringer by 
an array of line stiffeners. The stringer was divided 
into strips, and each strip was modeled by a line stiff- 
ener attached to the sheet at discrete points with no 
coupling between adjacent line elements. This anal- 
ysis did not model debond of the adhesive; thus, load 
predictions were unrealistically high in the line stiff- 
ener nearest the crack tip. 

In reference 7, the solution to this problem as- 
suming a linear elastic adhesive with an intact bond 
is presented. The parameter with the greatest in- 
fluence on the stress-intensity factors was found to 
be the distance from the crack tip to the edge of 
the stringer. Unless the crack tip was very close to 
or under the stringer, the stress-intensity factor was 
approximately that of an unstiffened sheet. How- 
ever, as the crack propagated beneath the stringer, 
the stress-intensity factor decreased significantly. In- 
creasing the stringer stiffness or the adhesive stiff- 
ness also resulted in a decrease in the stress-intensity 
factor. 

In reference 8, the solution to this problem assum- 
ing a nonlinear elastic adhesive with an intact bond 
is presented. When the adhesive was modeled with a 
nonlinear stress-strain curve, the peak stresses in the 
adhesive were considerably reduced in comparison to 
the solution for the linear elastic adhesive. This re- 
duction resulted in increases in the stress-intensity 
factors due to the nonlinearity of the adhesive. 


Experimental work has shown that as the crack 
tip approaches the stringer, debonding may start and 
propagate through the adhesive (ref. 5). Addition- 
ally, the adhesive may exhibit regions of nonlinear 
behavior. A need exists for an analysis that can be 
used to evaluate the effect of adhesive nonlinearity 
and a growing debond area on crack propagation in 
adhesively bonded structures. The purpose of the 
current work is to develop such an analysis and use 
it to predict the extent of stringer debonding and the 
effects of the debonding on the stress-intensity fac- 
tors for cracks in a composite-material sheet /stringer 
configuration. Trends are developed for the effects of 
applied stress level, location of the crack tip relative 
to the stringer edge, and linear versus nonlinear ad- 
hesive behavior. 

Nomenclature 


a 

half-length of crack, m 

+a 

right crack tip, m 

—a 

left crack tip, m 

b 

distance from edge of stringer 
to center of crack, m 

D 

domain of integration 

d X , dy 

constants in piecewise linear 
approximation of adhesive 
stress-strain curve 

Ex, Ey 

Young’s moduli in x- and in- 
directions, respectively, Pa 

/(Hi /()2 

stress functions, m 3 /N 

Gxy 

orthotropic shear modulus, Pa 

G 3 

shear modulus of adhesive 
layer, Pa 

hj 

thickness of layer j, m 

ki 

mode-/ component of stress- 
intensity factor, Pa^/m 

L/D 

linear elastic adhesive with 
adhesive debonding 

L/I 

linear elastic adhesive with 
intact bond 

N/D 

nonlinear elastic adhesive with 
adhesive debonding 

N/I 

nonlinear elastic adhesive with 
intact bond 

m x ,m y 

reciprocals of local slope in x- 
and y-directions, respectively, 
1/Pa 
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N number of cells or number of 

collocation points in domain 

R x , Ry major and minor axes of 

debonded area, m 

r x , r y coordinate axes of debonded 

area, m 

Sjk complex kernels in integral 

equations (j, k = 1,2), 
l/(Pa-m) 

Wi,W2 complex kernels in stress- 

intensity-factor equations 

Xj,Yj body forces acting on layer 

3 (3 = 1 . 2 ), N/m 3 


uniform pressure ao on the crack faces with the stress 
state at infinity equal to zero, as shown in figure 1. 
As explained in the appendix, the present solution for 
the uniformly stressed crack face is a close approxi- 
mation for the corresponding problem with a remote 
applied stress. 

Integral Equations for Adhesive Shear 
Stresses 

By using Green’s functions and the complex vari- 
able theory of orthotropic elasticity developed by 
Lekhnitskii, a set of integral equations describing the 
problem was formulated (refs. 7 , 8, and 9 ). For elastic 
nonlinear adhesive behavior, the problem in figure 1 
was reduced to the solution of the integral equations 
which follow. 


x, y Cartesian coordinates, m 

XQ,yQ coordinates of concentrated 

load point, m 

A Xj, Ayj incremental distances 

(. i = 1,N), m 

7 shear strain 

u X y, Vyx Poisson’s ratios 

ao crack-face pressure, Pa 

r shear stress, Pa 

Tj , 7 ,• arbitrary point on adhesive 

stress-strain curve 

tr rupture strength of adhesive in 

shear, Pa 

Tx, T y adhesive shear stresses, Pa 

Analytical Procedures 

Formulation of Problem 

Consider the configuration shown in figure 1 . It 
is assumed that the debond penetration under the 
stringer is small compared with the width of the 
stringer. Based on this assumption, the stringer is 
modeled as a semi-infinite sheet bonded to an infi- 
nite sheet with a finite, through-the-thickness crack 
in the infinite sheet. Hereinafter, the semi-infinite 
sheet is referred to as a stringer. The stringer and 
the infinite sheet are bonded together by an adhe- 
sive layer of constant thickness /13. It is assumed 
that the adherends are in a state of plane stress and 
the adhesive is in pure shear. Both the infinite sheet 
and the stringer are made of a fiber-reinforced plas- 
tic which is treated as a homogeneous, orthotropic, 
linearly elastic medium. The model is loaded by a 


h3m x T x {x,y) + j j [Sn(x,y, x 0 ,y 0 )r x (x 0 ,y 0 ) 

D 

+ {x, y, x 0 , t/o) T y (x 0 , 3/0)] dx 0 dy 0 

= <rofoi(x,y) - h 3 d x 

► 

h 3 m y Ty{x,y) + j j [S21 {x,y,xo,y 0 )T x (x 0 ,yo) 

D 

+ S22 (x, y, xq, yo) T y (xo, yo)] dx 0 dyo 
= <xofo2{x,y) - hzdy 

( 1 ) 

Here t x and T y are the unknown adhesive shear 
stresses in the x- and y-directions, respectively, and 
the kernels Sj^ ( j , k = 1,2) and the functions /01 and 
/02 are known; m x and m y are the reciprocals of the 
local slope of the adhesive shear stress-strain curve, 
and d x and d y are constant terms in the piecewise 
linear approximation of the adhesive stress-strain 
curve; D is the region of the adhesive bond. (The 
appendix details the derivation of equations ( 1 ).) For 
the linear adhesive case, equations (1) are simplified 
as follows: 

m x -m y = I/G3 


d x — d y — 0 

The solution of equations ( 1 ) produces the shear- 
stress distribution in the adhesive t x and r y . For the 
linear elastic representation of the adhesive, these 
equations need only be solved once. When the 
nonlinear stress-strain curve is used to model the 
adhesive, these equations are solved iteratively using 
a piecewise linear approximation of the adhesive 
stress-strain curve. Once the adhesive stresses are 
known, the stress-intensity factor at either crack tip 
can be found, and the debond area can be predicted. 
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Numerical Solution of Integral Equations 

Because of the complicated nature of equa- 
tions (1), a closed-form integration is difficult if not 
impossible to perform. Since the kernels Sjk (j, k = 
1, 2) contain only logarithmic singularities and are in- 
tegrate in the infinite domain D, equations (1) can 
be treated as Fredholm equations of the second kind. 
Consequently, the system of equations can be solved 
with standard numerical techniques. This is done by 
dividing the domain into small cells, assuming the 
unknown functions t x and T y to be uniform in each 
cell, and then using a numerical scheme to evaluate 
the equations. 

When the integrals in equations (1) are replaced 
by summations, the following system of equations is 
obtained: 

N 

him x T x (xj,yj) + [ S ll ( x j,yj, x On,yo n ) 

n=l 

x Tx (zo„ , 2/0„ ) + Sl2 ( x j> 2/j , x 0n i 2/On ) 

X Ty (zo„, 2 / 0 „)] Az„A y n 
= (Tofoi {xj, Vj) - h 3 d x ( j = 1, N) 

• ( 2 ) 

N 

hsmyTy ( Xj,yj ) + ^ [S 2 l {xj, yj, x 0n , yo n ) 

n= 1 

X T X (xo n , y 0 „ ) + s 2 2 (xj, Vj, X 0n , y 0n ) 
x T y (x 0n ,y 0n )\ Ax n A y n 
= Pof02 (Xj, Vj) - h 3 dy ( j = 1, N) 

where N is the number of collocation points or the 
number of cells in the domain. 

The outer boundary of the domain D theoreti- 
cally goes to infinity. However, to carry out the nu- 
merical analysis, the size of D must be restricted. A 
convergence study was conducted to determine the 
extent and refinement necessary so that the crit- 
ical quantities of interest, the stress-intensity fac- 
tors, were not appreciably affected by the restric- 
tion. Thus, for the linear elastic adhesive solution 
(ref. 7), the size of D was determined iteratively, 
starting with a small, coarse mesh and increasing the 
extent and refinement until no significant changes oc- 
curred in the stress-intensity factor. In reference 8, 
it was determined that the same mesh could be used 
with the nonlinear adhesive equations as was used 
with the linear adhesive. A typical mesh layout used 
in the analysis is shown in figure 2, where, because 
of the symmetry, only one-half of the integration' 
domain is shown. 


Determination of Stress-Intensity Factors 

In the present problem, because of symmetries 
in the geometry and loading, the shear component 
of the stress- intensity factor is zero. The normal 
component is found by combining the effects of the 
crack-face pressure uq and the adhesive shear forces 
t x and Ty acting on the infinite sheet. The stress- 
intensity factor may be expressed in terms of the 
unknown shear forces as follows: 

ki (oo) = crQ\fa + 

D 

+ W 2 (x,y,x 0 ,y 0 )Ty (x 0 ,yo)} d xo d yo ( 3 ) 

where a 0 is replaced by +a for the right crack tip 
and by —a for the left crack tip. The terms W\ and 
W 2 are functions of the crack length, edge distance 
of the stringer, and the material properties of the 
infinite sheet. Complete details on the derivation of 
the equation for the stress-intensity factors are given 
in reference 7. 

Prediction of Debonding 

To predict the boundary of the debond area, 
a simple rupture criterion based on the adhesive 
shear stresses is introduced. This criterion, shown 
in equation (4), states that within the debond area, 
the state of combined shear stress in the adhesive 
is greater than or equal to r^, the rupture strength 
of the adhesive in shear, where tr is known. The 
adhesive failure criterion is written as follows: 



This failure criterion is based on the octahedral 
shearing-stress theory (also referred to as the Mises- 
Hencky, or simply the Von Mises, criterion) where 
the ultimate shear strength of the adhesive tr has 
replaced the yield strength found in the Von Mises 
equation. Since the only stresses assumed to be 
present in the adhesive layer are the shear stresses, 
only these terms appear in the failure criterion. 

Linear adhesive. To determine the boundary of 
the debond area for the linear adhesive, the adhe- 
sive shear stresses in each cell are first found as- 
suming no debond. At the given applied stress, ev- 
ery cell where equation (4) is true is assumed to be 
debonded. Using this new debond area, the adhesive 
shear stresses are determined again and compared 
with equation (4) to find any additional debonding. 
If there is any, the new debond area is used to recal- 
culate the adhesive stresses. This iteration process 


[Wi ( x,y,x 0 ,yo)Tx (zq.Vo) 
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continues until there is no new debond area at the 
current level of applied stress. For each consecutive 
level of applied stress, the debond area determined 
from the previous level is used to start the iteration. 

Nonlinear adhesive. The debond area for the 
nonlinear adhesive is found in a similar manner. 
First, the adhesive shear stresses are found assuming 
no debonding of the adhesive layer. In this case, 
though, equations (2) must be solved iteratively. 
Then, as with the linear adhesive, every cell where 
equation (4) is true is assumed to be debonded; 
the procedure follows as in the linear adhesive case, 
where equations (2) are solved iteratively for each 
new assumed debond area. 

Results 

For the results presented, the infinite sheet 
(layer 1) is modeled as a quasi-isotropic graphite/ 
epoxy laminate and the stringer (layer 2) as a uni- 
directional graphite/epoxy laminate. The material 
properties (ref. 10) are given in table I. For the elas- 
tic linear adhesive, a shear modulus of 560 MPa was 
used. The rupture strength of the adhesive in shear 
tji is taken as 37.9 MPa. These values are typical 
of the AF-127 adhesive (ref. 5). The adhesive stress- 
strain curves used to model the adhesive behavior are 
shown in figure 3. The thickness of the base sheet hi 
is 2.0 mm, the thickness of the stringer h 2 is 1.0 mm, 
and the adhesive thickness h 3 is 0.01 mm. 

Effect of Stringer Debond Size and Shape on 
Stress-Intensity Factors 

The effect of adhesive debonding on the stress- 
intensity factors k\ is investigated by assuming var- 
ious debond sizes and shapes. Figure 4 shows the 
configuration of the assumed debond; R x and Ry 
are the half-lengths along the axes of the elliptical 
debond area in the x- and y-directions, respectively. 
The values of R x and Ry are varied to produce dif- 
ferent debond sizes. A value of zero for either Rx 
or Ry indicates no debonding of the adhesive (i.e., a 
zero debond area). Although the crack as shown in 
figure 4 is not beneath the stringer, for most of the 
results presented, the crack is partially beneath the 
stringer. In the following figures, only results for the 
crack tip nearest the stringer are shown, because the 
size of the debond area has almost no effect on the 
stress-intensity factors for the crack tip farthest from 
the stringer. 

Figure 5 shows a series of curves in which the 
value of R x is varied and Ry is held constant for 
a/b — 5.0. The location of the crack tip in relation 
to the edge of the stringer is noted in the figure. For 


all the values of Ry considered, the results show that 
the stringer must be debonded nearly to the end of 
the crack before the value of k\ begins to increase. 
Once the debond is extended past the crack tip, the 
maximum increase in k\ is quickly obtained, and 
further debond growth has little effect on h\. For 
increasing values of Ry, the height of the debond, 
the limiting value approaches 1.0, the solution for 
the unstiffened sheet, as expected. 

In figure 6, the value of Ry is varied while R x is 
held constant for a/b = 2.0. The debond width was 
chosen so that a/R x = a/b = 2.0; this locates the end 
of the debond area at the crack tip. As indicated 
in figure 6, in this example two debond shapes are 
considered: a rectangular area and an elliptic area. 
The change in the crack-tip stress-intensity factor 
as the debond areas increase is shown in the figure. 
The rectangular debond results in larger increases in 
k\ than does the elliptic debond. This is expected, 
since, for the same values of R x and Ry, a rectangle 
has more area than an ellipse. Thus, with the larger 
debond area, the stringer is less effective in reducing 
fci- However, as the value of Ry increases (i.e., as 
the debond grows along the edge of the stringer), 
differences between a rectangular debond and an 
elliptic debond decrease. As Ry approaches infinity, 
the solution for k\ should approach the solution for 
b = a, as indicated in figure 6. 

Figure 7 shows the effect of assuming that the 
debond has extended completely across the width 
of the stringer (i.e., R x = 00). The stress-intensity 
factor is plotted at various levels of Ry for increasing 
values of a/b. For values of a/b less than 1.0, the 
debond height has little effect on the crack-tip stress- 
intensity factor. However, for a/b > 1.0, that is, 
when the crack tip is beneath the stringer, the effect 
of the stringer is progressively reduced as Ry becomes 
larger. As Ry becomes very large, the values of k\ 
approach those of the unstiffened sheet. 

Predicted Stringer Debonding 

Linear adhesive. Figures 8, 9, and 10 show the 
predicted debond areas at increasing levels of applied 
stress for a/b = 0.95, 1.05, and 2.0, respectively. 
These values of a/b were chosen to illustrate three 
particular cases of interest: (1) the crack tip close to 
but not under the stringer; (2) the crack tip a small 
distance beneath the stringer; and (3) the crack tip 
a large distance beneath the stringer. 

For a/b — 0.95 (fig. 8), the right crack tip is 
located slightly beneath the edge of the stringer. In 
this case, debonding does not start until the applied 
stress reaches 70 MPa. With increasing applied stress 
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level, the debond grows significantly along the edge 
of the stringer but very little under the stringer. 

Figure 9 shows the predicted debond growth for 
a/b = 1.05. In this case, the crack tip is beneath 
the stringer, slightly past the edge of the stringer. 
The debonding begins at a much lower stress level 
for this case than for a/b = 0.95. For a/b — 1.05, 
the debond initiates at an applied stress of 15 MPa. 
Once the debond is initiated at the stringer edge, it 
grows first to the crack tip, and then the rest of the 
debond growth is primarily along the stringer edge. 

When the crack tip is far beneath the stringer 
(a/b = 2.0), as shown in figure 10, the debonding 
starts at an even lower applied stress than with the 
two previous cases. Here, the debond initiates at an 
applied stress of 7.5 MPa and grows diagonally until 
the debond reaches the crack tip. Once the debond 
reaches the crack tip, it continues to extend parallel 
to the edge of the stringer, as in the previous cases. 

Nonlinear adhesive. Using the stress-strain curve 
shown in figure 3, the debond growth for the nonlin- 
ear adhesive is predicted for a/b = 0.95, 1.05, and 
2.0. For the materials considered, the solution of the 
nonlinear equations typically required 8 to 10 itera- 
tions for convergence. 

Figure 11 shows the predicted debond area for 
a/b = 0.95, assuming a nonlinear elastic adhesive. 
The debond initiates at an applied stress of 150 MPa, 
which is much higher than the stress level required 
for debond initiation in the linear elastic adhesive 
case. Also, the debond does not start at y = 0 
(the crack plane), as in the linear adhesive case, but 
starts at the edge of the stringer some distance above 
the crack. Additionally, compared with the linear 
adhesive case, equal increments of stress produce less 
debond growth for the nonlinear adhesive. 

Figure 12 shows the predicted debond area for 
a/b = 1.05, assuming a nonlinear elastic adhesive 
layer. The debond initiates at an applied stress 
of 30 MPa, compared with a stress of 15 MPa re- 
quired to initiate the debond for the linear adhesive 
case. Comparison of figures 9 and 12 shows that the 
debond trends are similar for the linear and nonlinear 
adhesive models but that higher stresses are required 
to grow the debond for the nonlinear adhesive. 

Figure 13 shows the predicted debond area for the 
nonlinear adhesive for a/b = 2.0. Here, the debond 
initiates at an applied stress of 15 MPa, compared 
with 7.5 MPa required for debond initiation in the 
linear adhesive case. Again, comparison of figures 10 
and 13 shows that the debond does not grow as 
readily with the nonlinear adhesive as with the linear 
adhesive, although the trends are similar. 


Predicted Stress-Intensity Factors 

Figures 14 and 15 compare the normalized stress- 
intensity factors for four conditions: linear elastic ad- 
hesive with an intact bond (L/I), linear elastic adhe- 
sive including debonding (L/D), nonlinear elastic ad- 
hesive with an intact bond (N/I), and nonlinear elas- 
tic adhesive including debonding (N/D). Figure 14 
shows the normalized Aq for a/b = 1.05, and figure 15 
shows Aq for a/b = 2.0, at both the left, Aq(— a), and 
right, /q(+a), crack tips for increasing levels of ap- 
plied stress. For the intact conditions, in most cases, 
the applied stress levels used in these figures pro- 
duce adhesive stresses greater than r^, the adhesive 
rupture strength. However, for the intact cases, the 
adhesive was assumed not to rupture to allow for 
comparisons with the debond cases. Thus, for these 
data points, the adhesive stress-strain curves shown 
in figure 3 were extended past the point of adhesive 
failure in a straight line having the same slope as the 
final portion of the curves. 

As shown in figures 14 and 15, modeling adhe- 
sive debonding, nonlinearity, or the combination of 
the two, causes Aq to increase compared with the lin- 
ear intact case. For all conditions, except the linear 
intact case, normalized Aq increases with increasing 
level of applied stress. The effect of adhesive debond- 
ing and nonlinearity on the stress-intensity factors at 
the left crack tip is much less than at the right crack 
tip. 

Figures 14 and 15 show that, for the conditions 
examined, debonding of the linear adhesive produces 
the largest increase in Aq. Also, the greater the 
applied stress, the larger the increase in Aq. This 
is particularly evident for a/b = 2.0 (fig. 15). In this 
case, for an applied stress of 25 MPa, the debonding 
of the linear adhesive is so extensive that the effect 
of the stringer in reducing Aq is almost entirely 
eliminated. 

An interesting comparison can be made for the 
nonlinear adhesive with the intact and debond con- 
ditions from figures 14 and 15. (The data shown 
in figures 14 and 15 are listed in table II.) Includ- 
ing the nonlinearity of the adhesive with the intact 
bond condition causes the stress-intensity factors to 
increase relative to the linear intact condition. How- 
ever, including the nonlinearity of the adhesive with 
the debond condition causes the stress-intensity fac- 
tors to decrease relative to the linear debond con- 
dition. This agrees with the conclusion reached in 
reference 5, which was that modeling a nonlinear ad- 
hesive causes the stress-intensity factor to decrease. 
(The analysis of ref. 5 included debonding of the ad- 
hesive in the model.) This also supports the con- 
clusion of references 8 and 11, which showed that 


6 



modeling a nonlinear adhesive results in increases in 
the stress-intensity factors. (Neither refs. 8 nor 11 
considered debonding of the adhesive layer.) 

Conclusions 

This report presents the formulation of an anal- 
ysis for a cracked orthotropic sheet reinforced with 
an adhesively bonded orthotropic semi-infinite sheet. 
This configuration was assumed to represent a 
bonded stringer panel when the crack is located close 
to the edge of the stringer, and any debonding of the 
adhesive layer is small compared with the width of 
the stringer. The effect of debond growth and ad- 
hesive nonlinearity on the crack-tip stress-intensity 
factors was examined using the analysis. Based on 
this study, the following conclusions were drawn: 

1. For the linear elastic adhesive, debonding of the 
stringer near the crack caused the stress-intensity 
factors to increase relative to the intact bond 
condition, and, when the debond dimensions ex- 
ceeded the crack half-length, the stress-intensity 
factors approached those of the unstiffened sheet. 

2. When the adhesive was modeled as a linear elastic 
material, the analysis indicated that (a) when the 
crack is not beneath the stringer, the majority 
of the debond growth is along the edge of the 
stringer; (b) when the crack tip is beneath the 


stringer, the debond first grows to the crack tip 
and then along the edge of the stringer; and 
(c) the stress level required to start the debond 
growth in the adhesive layer decreases as the crack 
tip moves towards and beneath the stringer. 

3. When adhesive nonlinearity is modeled in the 
analysis, the applied stress levels required to ini- 
tiate debonding increase significantly compared 
with the stress levels required for the linear elas- 
tic adhesive, and the predicted debond does not 
extend as readily as for the linear adhesive. 

4. For both the linear and nonlinear representations 
of the adhesive layer, the stress-intensity factors 
increased when debonding was allowed to occur 
according to the assumed adhesive failure crite- 
rion. Including debonding with the linear elastic 
adhesive caused the stress-intensity factors to in- 
crease relative to the linear intact condition, but 
including debonding with the nonlinear adhesive 
caused the stress-intensity factors to decrease rel- 
ative to the nonlinear-intact condition. That is, 
when debonding was included, the stringer was 
more effective with the nonlinear adhesive. 

NASA Langley Research Center 

Hampton, VA 23665-5225 

October 9, 1985 
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Appendix 

Formulation of Integral Equations for 
Bonded System 

This appendix is a brief summary of the formula- 
tion of the integral equations for the solution to the 
problem of the semi-infinite orthotropic sheet bonded 
to a cracked orthotropic sheet. The complete details 
of this formulation can be found in references 7 and 9. 

Consider the stringer configuration shown in fig- 
ure 1. (Hereinafter, the semi-infinite sheet will be 
referred to as a stringer.) The stringer and the sheet 
are bonded together by an adhesive layer of constant 
thickness h 3. 

The model is loaded by a uniform pressure op on 
the crack faces with the stress state at infinity equal 
to zero as shown in figure Al. If the two layers have 
equal minor in-plane Poisson’s ratios u yx , the stress- 
intensity-factor solution formulated here for the uni- 
formly stressed crack face is identical to the problem 
of the remotely loaded stringer panel. On the other 

hand, if i>y~x 7^ Vyx, as is the case here, the displace- 
ments in the x-direction will differ for the two loading 
cases. Therefore, the two solutions will not be equiv- 
alent. However, it is shown in reference 9 that the 
differences in the two solutions are not large. For the 
graphite-epoxy laminates considered, the error intro- 
duced by this approximation is less than 5 percent. 
Thus, the present solution for a uniformly stressed 
crack face is a close approximation to the solution 
for the corresponding case with remote loading. 

For the problem shown in figure Al, the integral 
equations are formulated with the following assump- 
tions: 

1. The sheet (layer 1) and the stringer (layer 2) 
are homogeneous and linearly elastic; the adhesive 
(layer 3) is homogeneous, isotropic, and nonlinearly 
elastic. 

2. The sheet and the stringer are dissimilar, or- 
thotropic materials with principal directions of or- 
thotropy being oriented parallel and perpendicular 
to the crack in layer 1. 

3. The thickness of the sheet hi and the thickness 
of the stringer /12 are small compared with the in- 
plane dimensions, so that both layers are considered 
to be in a state of plane stress. 

4. The surface shear transmitted through the 
adhesive is assumed to act as a body force on the 
infinite sheet and the stringer. 

5. The thickness of the adhesive /13 is small com- 
pared with the thicknesses of the sheet and stringer; 
thus, the adhesive layer is treated as a shear spring 
rather than as an elastic continuum. 


Assumption 5 leads to the following continuity of 
displacement equations: 


h$ 

ui-u 2 = 77 -t x 


ho 

vi-v 2 = ^ rTy (Al) 


where u\, rq and U2, V2 are the x- and y-components 
of the in-plane displacement vectors in layers 1 and 
2, t x and r y are the components of the adhesive shear 
stress, and /13 and G3 are the adhesive thickness and 
shear modulus, respectively. From assumption 3, the 
two sets of body forces (force per unit volume) that 
act on layers 1 and 2 (see fig. Al) can be written as 
follows: 






(A2) 

Figure Al shows how the problem is broken into 
its component parts. The displacements in the sheet, 
shown in part B of figure Al, and in the stringer, 
shown in part C of figure Al, are determined individ- 
ually. The complete expressions for these displace- 
ments are given in reference 7. Equations (A2) are 
used to relate the forces X \ , Y\ , X2 , Y2 to the adhe- 
sive shear stresses t x ,t v . The displacements for each 
layer are substituted into equations (Al). Then, fol- 
lowing the procedure used in reference 11, a relation- 
ship between the adhesive shear stress and strain is 
obtained by linear interpolation. If r is the unknown 
stress on the stress-strain curve (fig. A2) and Tj-i 
and Tj are the adjacent points given on the stress- 
strain curve, where the corresponding shear strains 
are lj—i < 7 < 7y, then the relationship between 7 
and r can be written as follows: 


7 = 7j-l +{t~ Tj- 1 ) 


7 j ~ Tj-i 
T j ~ T j~ l 


7 - mr + d 


where 


Ij ~ lj-l 
m — — - — 

T j ~ T J~ l 

d = 7 j-i - rriTj-i 


and 


7 = r/G 3 
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Thus, for the nonlinear elastic adhesive, the prob- 
lem reduces to the solution of the following integral 
equations: 

h3m x T x {x,y) + J J [Sn{x,y,x 0 ,y 0 )T x {x 0 ,y 0 ) j 
D 

+ S12 (x, y, x 0 , ?/o ) Ty (xo, 2/0)] dxo dy 0 
= <rofoi{x,y) - h 3 d x 

h^m tJ Ty(x 1 y) + J j [S2i(x,y,x 0 ,y 0 )T x {xQ,y Q ) 


+ S22 y ' x 0 ’ Vo) T v ( x 0, yo)] dx 0 dy 0 
= <rofo2{x, y) - hzdy 


(A3) 


where D is the bonded region, and where the kernels 
Sjktiik = 1 > 2 ) and the functions /01 and /02 are 
known. (See ref. 7 for details.) The kernels Sji : are 
related to the distributed body force loadings shown 
in parts B and C of figure Al and are functions 
of the material properties of each layer, the half- 
crack length a, and the stringer edge distance b. 
The functions /qi and /02 are related to the uniform 
crack-face pressure a 0 shown in part B of figure Al 
and are functions of the material properties of layer 1 
and the half-length of the crack a. The terms m x and 
m y represent the reciprocal slopes of the adhesive 
stress-strain curve (defined previously) in the x- and 
y-directions, respectively, and d x and d y are constant 
terms (defined previously) in the piecewise linear 
approximation of the adhesive stress-strain curve. 
Complete details on the derivation of these equations 
can be found in reference 8 or 9. 
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A. Total system B. Sheet C. Stringer 

(layer 1) (layer 2) 



Figure Al. Superposition model of problem. 
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TABLE I. MATERIAL PROPERTIES 


Layer 


Ey, 

GPa 

Vxy 

G X y i 

GPa 

T R, 

MPa 

Material system 

1 

51.40 

51.40 

0.30650 

19.67 


[0/±45/90] s graphite-epoxy 

2 

10.86 

129.40 

.02617 

5.70 


[0] graphite-epoxy 

3 




.56 

37.9 

AF-127 adhesive 


TABLE II. PREDICTED STRESS-INTENSITY FACTORS 




Normalized stress-intensity factor at crack tip nearest stringer 1 for — 

a/b 

^0) 

MPa 

L/I 

L/D 

N/I 

N/D 

1.05 





0.7366 



.6679 


.7383 

.7725 


50.0 

.6679 

.8552 

.7522 

.7872 

2.00 

15.0 

0.2081 


0.2105 



20.0 


■ 

.2142 

■ 


25.0 

.2081 

.9163 

.2191 





Normalized stress-intensity factor at crack tip away from stringer 2 for — 

a/b 

<70> 

MPa 

L/I 

L/D 

N/I 

N/D 

1.05 

30.0 






40.0 

.9674 

.9697 

.9683 

.9687 


50.0 

.9674 

.9707 

.9685 

.9690 

2.00 

15.0 


0.8797 




20.0 

| | 

.8958 




25.0 

.8590 

.9965 


.8739 


1 fci(+a) 
cto s/'a ' 

2 fci(-a) 

ctq y/a 
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Figure 3. Shear stress-strain curve for AF-127 adhesive (ref. 5). 
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Figure 4. Configuration of debond area. 
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Figure 5. Normalized stress-intensity factor as a function of debond width and height, a/b = 5.0. 
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Figure 6. Normalized stress-intensity factor as a function of debond aspect ratio, a/b = 2.0; a/R x = 2.0. 
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Figure 9. Predicted debond growth in linear adhesive for a/b = 1.05. 






Figure 10. Predicted debond growth in linear adhesive for a/b = 2.0. 
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Figure 13. Predicted debond growth in nonlinear adhesive for a/b — 2.0. 
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